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ABSTRACT 

We report the discovery of periodic intensity dips in the X-ray source XMMU 
J004314.1-F410724, in the globular cluster Bol58 in M31. The X-ray flux was modu- 
lated by ~ 83% at a period of 2.78 hr (10017 s) in an XMM-Newton observation taken 
2002 Jan 6-7. The X-ray intensity dips show no energy dependence. We detected weaker 
dips with the same period in observations taken 2000 June 25 {XMM-Newton) and 1991 
June 26 {ROSAT/PSPC). The amplitude of the modulation has been found to be anti- 
correlated with source X-ray flux: it becomes lower when the source intensity rises. The 
energy spectrum of Bol58 was stable from observation to observation, with a characteris- 
tic cutoff at ~ 4 — 6 keV. The photo-electric absorption was consistent with the Galactic 
foreground value. No significant spectral changes were seen in the course of the dips. 
If the 2.78 hr cycle is the binary period of Bo 158 the system is highly compact, with a 
binary separation of ~ 10^^ cm. The association of the source with a globular cluster, 
together with spectral parameters consistent with Galactic neutron star sources, suggests 
that X-rays are emitted by an accreting neutron star. The properties of Bo 158 are some- 
what reminiscent of the Galactic X-ray sources exhibiting a dip-like modulations. We 
discuss two possible mechanisms explaining the energy-independent modulation observed 
in Bo 158: i) the obscuration of the central source by highly ionized material that scatters 
X-rays out of the line of sight; ii) partial covering of an extended source by an opaque 
absorber which occults varying fractions of the source. 



Subject headings: galaxies: individual (M31) — stars: individual (Bol58) — galaxies: 
star clusters — X-rays: galaxies — X-rays: stars 



1. INTRODUCTION 

The X-ray source XMMU J004314. 1+410724 
was discovered in M31 by the Einstein observa- 
tory (source #85 in Trinchieri & Fabbiano 1991) 
and was detected in subsequent observations with 
i?0S'.4r(Primini et al. 1993, Supper et al. 2001), 
XMM-Newton (Shirey et al. 2001) and Chandra 



(DiStefano et al. 2002). Based on the Chan- 
dra aspect solution, which is currently limited by 
systematics to ~ 0.6" accuracy, the source loca- 
tion is a = 00''43'"14.42^ 6 = 41°07'26.3" (2000 
equinox) (DiStefano et al. 2002) (Fig. 1). This 
position of the source is consistent with the opti- 
cally identified globular cluster candidate Bo 158 
(source #158 in Table IV of Battistini et al. 1987) 



^Also: Space Research Institute, Russian Academy of Sciences, Profsoyuznaya 84/32, Moscow 117810, Russia 
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(we will use the designation "Bo 158" as a source 
name throughout this Letter). 

The previous observations were too insensitive 
for detailed study of individual sources in M31. 
The large collecting area and bandpass of XMM- 
Newton allow us to study the short-term variabil- 
ity and spectral properties of these sources. In 
this Letter wc report the discovery of periodic X- 
ray modulation of the lightcurve of Bo 158 and 
report its spectra. 

2. OBSERVATIONS AND DATA ANALYSIS 

In the following analysis we use data from three 
XMM-Newton observations of the bulge of M31 
(Table 1, Figure 1). The first observation of 
the central part of M31 was performed on June 
25, 2000 as a part of the Performance Verifica- 
tion Program (PI: M.G. Watson) (Shirey et al. 
2001, Osborne et al. 2001). Two other observa- 
tions were performed on June 29, 2001 (Shirey 
2001) and on January 6, 2002 as a part of the 
Guaranteed Time Program (PI: K.O. Mason and 
M.G. Watson). Wc use data from three Eu- 
ropean Photon Imaging Camera (EPIC) instru- 
ments: two EPIC MOS detectors (Turner et al. 
2001) and the EPIC PN detector (Strueder et al. 
2001). In all observations the EPIC instruments 
were operated in the full window mode (30' diam- 
eter FOV) with medium (2000 June 25 and 2001 
June 29 observations) and thin (2002 January 6 
observation) optical blocking filters. 

During all three XMM-Newton observations, 
the X-ray source Bo 158 was offset by ~ 10' from 
the center of the field of view (Fig. 1). In spite 
of the significant degradation in the sensitivity of 
the EPIC cameras at high offset angles, the statis- 
tics (more than 3000 counts in each observation) 
were sufficient for a detailed spectral and timing 
analysis. 

We reduced EPIC data with the XMM-Newton 
Science Analysis System (SAS v 5.3)^. Wc per- 
formed standard screening of the EPIC data to 
exclude time intervals with high background lev- 
els. To generate hghtcurves and spectra of the 
source, wc used an extraction radius of ~ 60" and 
subtracted as background the spectrum of adja- 
cent source-free regions with subsequent normal- 
ization by a ratio of detector areas. We used data 
in the 0.3 — 10 keV energy band because of the 
^Sec http: / /xmm. vilspa.esa.es / user 



uncertainties in the calibration of the EPIC in- 
struments outside this range. All fluxes and lu- 
minosities presented below apply to this band. In 
the following analysis we assume a source distance 
of 760 kpc (van den Bergh 2000). 

We used standard XANADU/XRONOSv.5 ^ 
Fourier transform and epoch folding tasks to 
search for the periodic modulation of the source 
X-ray flux and determine its period. 

We used spectral response matrixes generated 
by SAS tasks. The energy spectra of the source 
were fitted to two analytic models using XSPEC 
v.ll (Arnaud 1996): an absorbed simple power 
law (powerlaw) and comptonization (comptt) 
models^ EPIC-PN, MOSl and M0S2 data were 
fitted simultaneously, but with independent nor- 
malizations. 



3. RESULTS 

We generated X-ray lightcurves of the source in 
the 0.3 — 10 keV energy band using data from the 
EPIC-PN, MOSl and M0S2 detectors, combining 
them in order to improve statistics. The resulting 
hghtcurve of Bo 158 for the 2002 Jan. 6 observa- 
tion (Fig. 2 a) shows a pattern of recurrent dips 
in the X-ray intensity with a period of ~ 10000 
s. The dips are broad, lasting for ~ 30% of the 
10000-s cycle. Their FWHM varies between 2200 
and 3000 s, and their depths range from ~ 80 to 
~ 100% of the out-of-dip flux. The recurrence pe- 
riod is 10017 ± 50 s. Figure 26, c shows the power 
density spectrum of the source from 2 x 10~^ to 
3 X 10^^ Hz, and our determination of the pe- 
riod with an epoch folding analysis. Figure 2d 
shows the combined EPIC (PN+M0S1+M0S2) 
light curve for the Jan. 6 observation folded on 
the 10017 s best period of the modulation (phase 
0.0 is set arbitrarily at TJD=12279.0). In order 
to investigate the energy dependence of the X-ray 
modulation during 2002 Jan. 6 observation, we 
constructed lightcurves in the 0.3 — 2.0, 2.0 — 5.0, 
and 5.0 — 10.0 keV bands. As it is clearly seen 
in Figure 3 a, there is no significant energy depen- 
dence of the modulation. 

We searched for this modulation in earlier ob- 
servations of M31, including 2000 Jun. 25 and 
2001 Jun. 29 XMM-Newton observations of the 
center of M31, and ROSAT/PSPC observation 



'^http://heasarc.gsfc.nasa.gov/docs/xanadu/xronos/xronos.html 

^For model description see littp://lieasarc. gsfc.nasa.gov/docs/xanadu/xspec/xspec.litml and references therein 
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from 1991 Jul. 26. We produced combined EPIC- 
MO S lightcurves in the 0.3 — 10.0 keV energy 
range and a PSPC lightcurvc in the 0.2 — 2.0 
keV energy band. X-ray flux modulations with 
periods close to ~ 10000 s (~ 10500 s and 
~ 9700 s respectively) were marginally detected 
in the 2000 Jun. 25 and 1991 Jul. 26 observa- 
tions. Figure 3 shows the XMM-Newton/EPlC 
and ROSAT/PSPC light curves folded on the 
10017 s period determined above. The lightcurves 
of 2000 Jun. 25 and 1991 Jun. 26 observations 
show a ~ 30% and ~ 50% intensity drop during 
the dip (Fig. 3b, c), while for the 2001 Jun. 29 ob- 
servation, dips were not detected with a 2a upper 
hmit of 10%. 

In case of 2002 Jan. 6 and 2000 Jun. 25 observa- 
tions, where regular dips in the X-ray lightcurve of 
the source were clearly detected, we fit the folded 
0.3 — 10 keV lightcurves with a simple model, con- 
sisting of a constant plus a Gaussian with nega- 
tive normalization centered at the dip minimum: 
f(p) = C-Axexp[-(p-pdipy/2a^)] '^(thc best-fit 
model approximation is shown with dotted lines 
in Fig. 2d and Fig. 36). Using this approxi- 
mation, we obtain a 83 ± 5% amplitude of the 
dip with respect to the out-of-dip intensity and 
a = 0.10 ± 0.01 (in phase units of 10017 s cy- 
cle) for the 2002 Jan. 6 observation, and 30 ± 9% 
amplitude and a = 0.07 ± 0.03 for the 2000 Jun. 
25 observation. The values of a imply an average 
FWHM for the dips of order 2500 s. 

The averaged XMM/EFIC spectra of Bo 158 
was analyzed by fitting two different spectral 
models (Table 2). A model with a quasi- 
exponential cut-off (the Comptonization model in 
Table 2) at ~ 4.0 — 6.0 keV describes the energy 
spectra significantly better than a simple power 
law. This spectrum is somewhat reminiscent of 
Galactic neutron star systems with high luminos- 
ity (laria et al. 2001, DiSalvo et al. 2001). The 
spectrum is stable from observation to observa- 
tion, despite the change in X-ray luminosity from 
~ 6.5 X 10^7 (2002 Jan. 6) to ~ 1.4 x lO^^ ergs 
s~^ (2001 Jun. 29), although the cut-off energy 
may increase at the higher luminosities. 

In order to investigate any changes in the en- 
ergy spectrum during the intensity dips, we di- 
vide the spectra from 2002 Jan. 6 observation 



into two states corresponding to the out-of-dip 
and dip intensity (Fig. 4) (i.e. the fluxes higher 
and lower than 0.1 cnt/s in Fig. 2 a). A saturated 
comptonization model with electron temperature 
kTg ~ 1.7 keV and optical depth r ~ 22 gives 
a good approximation to both spectra (Fig. 4). 
The resulting 2(t upper limit on an increase in 
neutral absorbing column density A^h during the 
dips is 2 X 10^^ (assuming a single spectral com- 
ponent). However, if the dips are caused by elec- 
tron scattering in a cloud along the line of sight, 
the electron column density must be ~ 2.7 x 10^^ 
cm-2 (83% dipping) and ~ 5.4 x 10^=^ cm-^ (30% 
dipping) for 2002 Jan. 6 and 2000 Jun. 25 obser- 
vations respectively. 

It should be mentioned that the Comptoniza- 
tion model which flts the spectrum during the in- 
tensity dips does not provide a unique determina- 
tion of the spectral form. Other complex models 
(e.g. two-component models with different low- 
energy absorption column densities) also give ac- 
ceptable flts. The low statistical signiflcancc of 
the data does not allow us to discriminate be- 
tween these models. 

There is an anti-correlation between the aver- 
age amplitude of the modulation and the source 
intensity on a long time scale. The values of mod- 
ulation fraction were 83% in 2002 Jan. 6 observa- 
tion, 30% in 2000 Jun. 25 observation and < 10% 
in 2001 Jun. 29 observation, corresponding to the 
source luminosities of 6.5 x 10'^'', 9.9 x 10^^ and 
1.4 X 10^^ ergs s~^ respectively. In addition, a de- 
pendence of the dip strength upon the luminosity 
on a time scale of hours was marginally detected 
during the 2000 Jun. 25 observation. 

4. DISCUSSION 

We have discovered a 2. 78- hour periodic dip-like 
modulation with variable amplitude in the X-ray 
flux of the globular cluster candidate source Bo 
158 in M31. The amplitude of the modulation is 
anticorrelated with source intensity: it becomes 
lower with increasing X-ray flux. The most in- 
teresting feature of this modulation is its lack of 
energy dependence. 

If the 2.78-hour modulation represents the bi- 
nary period, Porb, the binary separation is a ~ 
7 X 10^°Mx^^(l + qy/^ cm ^ The association of 



^where p is a phase of 10017-s cycle, C is a source out-of-dip intensity level, A, pcUp and 2.35a are normalization, 
centroid phase and FWHM of Gaussian respectively 

^where Mx is a mass of the compact object in solar units and g is a mass ratio of the secondary star and a compact 
object 
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the source with a globular cluster, and spectral 
parameters consistent with Galactic neutron star 
sources, suggest that the compact object in Bol58 
is probably a neutron star in a low-mass binary. 
This implies a highly compact binary with sepa- 
ration a ^ 10^^ cm. 

The properties of Bo 158 are somewhat remi- 
niscent of Galactic X-ray sources that exhibit dip- 
like modulations (e.g. Gygnus X-3 and some "dip- 
ping" sources). 

Although there are similarities between the 
light curve of Bo 158 and the Galactic binary 
Gygnus X-3, it is unlikely that these systems have 
a common nature. The broad minimum in the 
4.8-hour cycle of Gygnus X-3 is explained as vari- 
able scattering in the dense photoionized wind of 
a Wolf-Rayet companion (Paerels et al. 2000). 
Moreover, in contrast to Bo 158, the amplitude 
of X-ray modulation in Gygnus X-3 is extremely 
stable in spite of drastic changes in luminosity. 

The duration of the intensity dips in phase is 
similar to those seen in some Galactic LMXB 
dipping sources such as X1755-338 (Mason ct al. 
1985), XB1254-690 (Gourvoisier et al. 1986) and 
XB 1916-053 (Yoshida et al. 1995), although Bo 
158 is much more luminous. The lack of energy 
dependence of the modulation is a common fea- 
ture of a number of dipping sources. The dips in 
these systems are believed to be caused by absorp- 
tion and scattering in an obscuring medium, i.e. 
a bulge or thickened region of the accretion disk, 
or inhomogeneities in the mass transfer stream 
(White & Swank 1982, Lubow & Shu 1976, Prank, 
King & Lasota 1987). 

If the dips in Bo 158 are caused by an obscur- 
ing structure, at least two explanations could be 
proposed to explain the energy independence of 
the intensity dips: i) the obscuring material is 
highly ionized (Mason et al. 1985); ii) the dips are 
caused by partial covering of an extended source 
(Ghurch et al. 1997). 

(i) If the obscuring structure is located close to 
the X-ray emission region, the obscuring medium 
could be strongly ionized. The compactness and 
high X-ray luminosity of Bo 158 could imply 
a high ionization in the region responsible for 
the dips. The obscuring matter in the Galactic 
dipping sources, with their smaller luminosities 
(10^^ — 10^^ ergs s~^), is less strongly ionized. For 



the material causing the dips to be completely 
ionized, the ionization parameter ^ = Lx/nR'^ 
(where Lx is the central source luminosity, n is 
the gas density of the cloud, and R is the dis- 
tance from the central source to the obscuring 
medium) (Hatchett, Buff & McGray 1976) should 
be larger than ~ 1000 ergs cm s~^. For the 2002 
Jan. 6 observation, an average 83% intensity re- 
duction requires a column, A'^e, of ~ 2.7 x 10^'' 
cm~^. Assuming an out-of-dip luminosity of Bo 
158 to be 6.5 x 10^'' ergs s~^, and following Ma- 
son et al. (1985), we conclude that all material 
closer than ~ 2.4 x 10^°e cm ^ will be highly ion- 
ized. Another estimate of the ionization stage 
of the obscuring medium is based on the dura- 
tion of the intensity dips (Remillard & Ganizares 
1984). For t^ip ~ 2500 s, assuming that the ab- 
sorbing medium is roughly spherical and fixed in 
the frame of the binary, ^ = 27rLxtdip/(A^-Porb-R)- 
To obtain ^ ^ 1000, the ionized region must be 
closer than i? ~ 4 x 10^° cm (2002 Jan. 6 obser- 
vation) and ~ 2 X 10^^ cm (2000 Jun. 25 ob- 
servation). The absorbing medium in a 2.78-hour 
binary could fall well within these limits, and be 
highly ionized. 

(ii) For an opaque absorber to produce broad 
dips, the emission region must be extended 
and the eclipse partial. The source could be 
a scattering corona-like structure, and/or thick 
disk/outflow located close to the compact object, 
and the absorber a thickened region of accretion 
disk, or a mass transfer stream from the compan- 
ion star. 

The diminution of the modulation with increase 
of the X-ray flux could indicate a decrease of the 
obscured fraction of the emission region, caused 
by increasing the size of the emitter or decreasing 
the size of the absorber. 

We would like to thank the referee for his/her 
helpful comments. This paper is based in part 
on observations obtained with XMM-Newton, an 
ESA science mission with instruments and con- 
tributions directly funded by ESA Member States 
and the USA (NASA). This research has made use 
of data obtained through the High Energy Astro- 
physics Science Archive Research Genter Online 
Service, provided by the NASA Goddard Space 
Fhght Genter. 



''where filling factor, e = {l/R), I is the scattering column length 
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Table 1 

XMM-NEWTON observations of M31 USED IN THIS ANALYSIS. 



Date, UT 


Tatart , 


UT 


Field 


Obs. ID 


RA (.72000)" 


Dec (,72000)° 


Exp.(MOS)'' 


Exp.(PN)'' 




(h;m 








(h:m:s) 


(d:m:s) 


(ks) 


(ks) 


25/06/2000 


10:44 


42 


M31 Core 


0112570401 


00:42:43.0 


41:15:46.1 


34.8 


31.0 


29/06/2001 


06:21 


38 


M31 Core 


0109270101 


00:42:43.0 


41:15:46.1 


32.6 


30.8 


06/01/2002 


18:07 


17 


M31 Core 


0112570101 


00:42:43.0 


41:15:46.1 


63.0 


61.0 



— coordinates of the center of the FOV 



— instrument exposure used in the analysis 



Table 2 

Best-fit model parameters of the energy spectra of Bo 158 during June 25, 2000, June 29, 
2001 and January 6, 2002 XMM - Newton observations of the central part of M31 (combined 
EPIC-PN, MOSl AND MOS2 data). Parameter errors correspond to la level. 

Parameter Observation Date 

Jun. 25, 2000 Jun. 29, 2001 Jan. 6, 2002 



Absorbed Power Law (powerlaw*wabs) 



Photon Index 


66+^-"^ 

^•""-0.03 




0.66 ± 0.03 




0.09 ± 0.02 


0.05 


-0.01 


Flux'^ 


1.738 ±0.035 


2.401 ±0.031 


1.082 ±0.019 


X'(d.o.f) 


302.3(262) 


414.6(355) 


287.4(267) 



Absorbed Comptonization Model (comptt*wabs) 





0.08+H^ 


0.08+H^ 


omtro'r 


kTd 


1 70+0.10 


1.87;°:°^ 


1.73 ±0.07 




19.7+li 


2o.il?:J 


21.6+12 


Nfj 




n 05+°-°^ 

'-'•^'3-0.05 




Flux'' 


1.442 ±0.029 


2.005 ± 0.026 


0.914 ±0.016 


X'(d.o.f) 


263.8(260) 


358.1(353) 


194.4(265) 



^ - equivalent absorbing hydrogen column density in units of 10 cm~ 

^ - model flux in the 0.3 — 10.0 keV energy range in units of lO"-*^^ erg cm."^ 

- temperature of soft photons, keV 

- electron temperature, keV 

^ - Thomson optical depth for the spherical geometry 
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Fig. 1. — X-ray image of the central part of M31, as it appears in the 63-ks exposure with EPIC-MOSl in 
the 0.3 - 10 keV energy range. The position of the X-ray source XMMU J004314. 1-1-410724 associated with 
globular cluster Bol58 (Battistini et al. 1987) is marked with an arrow. 
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Fig. 2. — (a) X-ray lightcurve of Bol58 during Jan. 6, 2002 XMM-Newton observation, obtained from 
combined data of EPIC-PN, MOSl and M0S2 cameras, 0.3 — 10 keV energy range, 500 s time resolution. 
(&) Power density spectrum of Bo 158. (c) Resulting distribution for the epoch folding analysis of the 
X-ray lightcurve. The best-fit value of the period (10017 s) is shown by a vertical dotted line, {d) Normalized 
lightcurve of the source, folded at the best-fit period of 10017 s. The analytic approximation is shown with a 
dotted line (see text). 
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Fig. 3.— X-ray lightcurves of Bol58 taken with XMM-Newton/EFIC and ROSAT/PSPC, folded on a period 
of 10017 s. (a) Folded X-ray lightcurves of Bol58 in the 0.3 — 2, 2 — 5 keV and 5 — 10 keV ranges during the 
January 6, 2002 XMM-Newton observation are shown with dotted, long-dashed and thin solid lines respectively 
(combined EPIC-MOS and FN data), (b) Folded X-ray lightcurve of Bol58 in the 0.3 — 10 keV energy range 
during the 2000 June 25 XMM observation (combined EPIC-MOS and FN data). The analytic approximation 
is shown with dotted line (see text), (c) Folded X-ray lightcurve of Bol58 in the 0.2 — 2.0 keV energy range 
during the 1991 June 26 ROSAT/PSPC observation. 
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Fig. 4. — The energy spectra of Bol58 from the non-dip (high flux) and dip (low flux) intervals of X-ray 
flux during the 2002 Jan. 6 observation. EPIC-PN data in the 0.3 — 10 keV energy range. The best-fit 
analytic models are shown with red and blue histograms, and are both saturated comptonization models with 
an electron temperature of ~ 1.7 keV and an optical depth of ~ 22. 



